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A range of monodomain nematic liquid crystal elastomers containing differing proportions of 
photo-isomerisable mesogenic moieties, which turn from a rod-like to a kinked shape upon ultraviolet 
(UV) irradiation, was studied. Depending on the proportion and positional role of the photo- 
sensitive groups in the crosslinked polymer network, different types and magnitudes of response were 
found. The principle consequence of such photo-isomerisation is the destabilisation of the nematic 
phase, whose order parameter depends on temperature in a near-critical fashion. Accordingly, the 
effect of UV-irradiation is dramatically enhanced near the critical temperature, with the associated 
reduction in the nematic order parameter manifesting as a change in the macroscopic shape of the 
elastomer samples, producing a large uniaxial contraction. Theoretical analysis of this phenomenon 
gives a good quantitative agreement with experiment. 

PACS numbers: 61.30.-v, 82.50.Hp, 61.41. +e, 82.30.Qt 



I. INTRODUCTION 

The unique behaviour of liquid crystal elastomers de- 
rives from the intimate relationship between the elastic 
nature of the polymer network, and the ordering of their 
mesogenic, liquid crystalline moieties. Many interest- 
ing and unusual properties of these materials have been 
identified, and are summarised in recent review articles 
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The property that is most relevant to this paper is the 
response of a liquid crystal elastomers macroscopic shape 
to temperature. A permanently, uniaxially aligned mon- 
odomain nematic liquid crystal elastomer network || will 
exhibit a spontaneous contraction along its director axis 
when heated towards its nematic-isotropic phase transi- 
tion temperature |t[ |8|. This phenomenon is due to the 
coupling between the average polymer chain anisotropy, 
and the nematic order parameter Q(T) ||. The large 
variation in magnitude of resulting spontaneous uniaxial 
deformation depends strongly on the extent of this cou- 
pling between the elastic and the liquid-crystal degrees 
of freedom in the material J^, . 

Accordingly, a mechanical action can be produced in 
response to a control signal which defines the underlying 
order parameter of the elastomer. Substantial macro- 
scopic movement and shape changes are brought about 
by subtle variations in microscopic ordering. There are 
many methods of achieving mechanical motion in re- 
sponse to a control signal - or actuation - for many 
different applications plf , ranging from precise piezo- 
electric positioning devices to artificial muscles jl2|, [l3| . 
Liquid crystal actuators predominantly rely on thermal 
|l(], Q and, to a much lesser extent, electromagnetic 
control signals. Recently, however, a less invasive 
method of manipulating the underlying nematic order 
parameter in an elastomer network has been reported for 
the first time 

When exposed to ultraviolet (UV) radiation, liquid 
crystal molecules containing photo-isomerisable groups, 
such as azobenzene, experience a reduction in their ne- 
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FIG. 1: Schematic depiction of nematic-isotropic phase trans- 
formation in a liquid crystal containing photo-isomerisable 
mesogenic molecules, which turn from a rod-like trans to a 
kinked cis conformation under the UV irradiation. 



matic order and can therefore be induced to undergo a 
phase transformation from the nematic to the isotropic 
state |l7| |l^]. This is due to the UV-induced cis- 
trans isomerisation of the azo- (N=N) bond, whereby 
the photo-sensitive mesogenic molecules change from a 
rod-like shape, to a strongly kinked shape upon irradia- 
tion at a resonant wavelength of 365nm [[l5[ [l?], Q . The 
rod-like shape serves to stabilise the liquid crystal phase, 
whereas the kinked shape has the opposite effect, acting 
as an impurity and destabilising the nematic phase, re- 
ducing its order parameter jnj, Fig. |l]. The dynamics of 
this process, and the reverse transacts reaction occur- 
ring on heating or irradiation at 465nm, are well studied 
in ordinary low- molar weight liquid crystals 0, p0[ . 

With this method of non-invasive control over the de- 
gree of nematic ordering, there is now a new mechanism 
for actuation in liquid crystal elastomers. Such mechan- 
ical photo-induced actuation in nematic elastomers was 
first reported by Finkelmann et al. [ p^5[ . The degree 
of wide- spread, far-reaching interest in this new photo- 
mechanical effect is highlighted by recent popular articles 
in the New Scientist and The Financial Times , which 
appear even before any in-depth research had been con- 
ducted, and suggest exotic applications such as surgical 
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micro-actuators. 

In this article, a range of such photo-sensitive nematic 
liquid crystal elastomers is studied. Before investigating 
their varied response to ultraviolet light, their thermal 
responses are first characterised in section [HI A , The dif- 
fering UV responses for the di fferent sample compositions 
are presented in section IIIB. Section^ then focuses on 
the detailed quantitative analysis of the UV-mechanical 
response, comparing th e resu lts with theoretical expres- 
sions derived in section IV B| . 



II. EXPERIMENTAL 

All starting materials and resultant aligned, mon- 
odomain nematic liquid crystal elastomers were prepared 
in the Cavendish Laboratory. The procedure for forming 
the side-chain polysiloxanes by the hydrosilation of the 
terminal vinyl group in the mesogenic moiety with the 
Si-H bond of the polysiloxane chain, as well as the two- 
step crosslinking technique with a uniaxial stress applied 
after the first stage of crosslinking to produce and freeze 
the monodomain nematic alignment, has been developed 
over the years by Finkelmann et al. |], ||, p2[ . A num- 
ber of minor modifications were made to the procedure, 
which are described in [10 . The azobenzene compounds 



were all synthesised according to standard literature tech- 
niques |23|, |24|, [25| . A summary of the different mesogenic 
moieties used is given in Fig. 

From the structure of the molecular moieties one can 
deduce that all the side-group compounds are mesogenic 
rod-like units. PBB and BMAzo have a four carbon 
atom spacer, which induces a parallel orientation of the 
rods to the siloxane backbone resulting in a prolate chain 
anisotropy in the nematic phase, with its principal radii 
of gyration i?y > Rj_ . The third rod-like side group, 
UBAzo, has a flexible spacer of eleven carbon atoms and 
is thus only expected to be weakly coupled to the silox- 
ane backbone. There are three distinct types of di- vinyl 
crosslinking agents: the flexible, no n- mesogenic DiUB, 
the rod-like DiUAzo and the main-chain nematic poly- 
mer DiMC, which contains a total of approximately 75 
rod-like monomer units between its reacting vinyl groups. 
The polymer backbone onto which all these compounds 
were grafted was polyhydromethylsiloxane (PHMS), con- 
taining approximately 60 SiH groups per chain, obtained 
from ACROS Chemicals. 



The polymer networks were all crosslinked via the same 
hydrosilation reaction, in the presence of a commercial 
platinum catalyst COD, obtained from Wacker Chemie. 
The crosslinking density was 10 mol% (molar percent, 
relative to the total number of SiH units) in all resultant 
elastomers. 

Phase sequences were established on a Perkin Elmer 
Pyris 1 differential scanning calorimeter, which corre- 
lated with the critical temperatures obtained by the ther- 



mal expansion measurements (section [III A ). In all re- 
sulting elastomeric materials a broad nematic liquid crys- 
talline phase below a certain critical temperature was ob- 
served, with optical microscopy and mechanical testing 
confirming that this was indeed the nematic phase. Some 
of the materials also possessed a second ordered phase at 
low temperatures, exhibiting high rigidity and little re- 
sponse to temperature or UV-irradiation, just above the 
glass transition. This additional phase was not a focus of 
the present work - all measurements were performed in 
the soft rubbery nematic phase of the elastomer samples. 

Elastomer mechanical expansion measurements were 
obtained by freely suspending the samples in an active, 
glass-fronted Gallenkamp oven and measuring their vari- 
ation in length with a travelling microscope. Cooling and 
heating rates were chosen such that thermal hysteresis ef- 
fects were negligible (around 0.5°C/min, except for the 
sample containing DiMC, which required a much slower 
rate). 

Samples (of approximate dimensions 30 x 5 x 0.4mm) 
were uniformly irradiated with a Merck VL-4.L UV lamp, 
providing a narrow band at 365nm with an output of 
approximately 4W, at a distance of approximately 15cm. 
The temperature was determined with a thermocouple 
placed directly behind the samples. 



III. COMPARISON OF MATERIAL 
PROPERTIES 



A wide range of elastomer compositions and topologies 
(summarised in the Table |) was investigated, in order to 
probe different aspects of nematic elastomer response to 
UV-irradiation. It was universally found that the main 
consequence of UV-irradiation was the cis-trans isomeri- 
sation of the azobenzene groups, with a very low prob- 
ability of side reactions. However, the differing compo- 
sitions and topologies of the samples resulted in a range 
of substantially different magnitudes of the effect. One 
key factor to consider is the overall concentration of iso- 
merisable azobenzene groups, which varied between 9 and 
100mol% in our materials. Another crucial factor is the 
topological role of the crosslinker (cf. Fig. ^]). The small, 
flexible crosslinking polymer, DiUB is deemed to have a 
minor effect on the overall mesogenic properties of the liq- 
uid crystal elastomer. This is quite unlike the other two 
crosslinkers - DiMC is a highly anisotropic main-chain ne- 
matic polymer containing ~75 rod-like monomer units, 
which can align itself and fold into hairpin defects in the 
nematic phase. DiUAzo contains a single, rod-like photo- 
isomerisable unit of its own; its effect on the mechanical 
response of the elastomer network to UV-irradiation was 
the subject of the pioneering study by Finkelmann et al. 

Pi. 
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FIG. 2: The mesogenic moieties used are abbreviated as follows: PBB - 4-pentylphenyl-4'-(4-buteneoxy)benzoate; 
BMAzo - [4-(4-buteneoxy)-4'-methyloxy]azobenzene; UBAzo - [4-(ll-undeceneoxy)-4'-butyloxy]azobenzene; DiUB - di- 
l,4-(ll-undeceneoxy)benzene; DiMC - a-4-[l-(4'-ll-undeceneoxybiphenyl)-2-phenyl]butyl-a;-(ll-undeceneoxy) poly-[l-(4- 
oxydecamethyleneoxy)-biphenyl-2-phenyl] butyl (~75 units long); DiUAzo - di-[4-(ll-undeceneoxy)]azobenzene. 



TABLE I: Compositions, and corresponding abbreviations, of the different photo-sensitive nematic elastomers studied in this 
work. The composition is presented in two forms: first, in mol% of corresponding moieties; second, the proportion of rod- like 
mesogenic units in the overall composition. The rod-like side groups are a combination of a "standard" nematic PBB and one 
of two kinds of the azo- containing groups (cf. Fig. ||), specified in the second column of the table. The crosslinking density in 
all cases is kept at a constant 10mol%. 
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Side j 


Composition (mol%) 
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10 
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9 
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A. Uniaxial thermal expansion 

Before their response to ultraviolet light was inves- 
tigated, the spontaneous uniaxial thermal expansion of 
each of the samples was first characterised (Fig. ||). This 
determines the extent of the coupling between the lo- 
cal nematic order, described by the scalar order parame- 
ter Q(T), and the average elastomer network anisotropy, 
measured as the spontaneous uniaxial extension A = 
L/Lq, where L is the current length of the sample, and 
Lo its length in the isotropic phase (the latter is always 
shorter for an elastomer with a prolate backbone config- 
uration, as is the case for all our materials). 

Empirically, one can relate the uniaxial thermal expan- 
sion of a liquid crystal elastomer, A, to the underlying 
nematic order parameter Q(T), which can be indepen- 
dently determined by either wide-angle X-ray scattering, 
or from the optical birefringence of the elastomer. It 
has been shown [pi p2l that the relationship between the 



macroscopic shape and the nematic order of liquid crystal 
elastomers can be very well described by a direct linear 
proportionality, A = 1 + aQ. The constant a is a sensi- 
tive function of the average network anisotropy, varying 
from a ~ 0.08 in weakly anisotropic elastomers to a ~ 3 
or greater in materials containing main-chain crosslinkers 
. Such a relationship between the spontaneous uniax- 
ial extension and the underlying nematic order is well- 
understood theoretically within the molecular model of 
equilibrium nematic polymer networks ||. 

However, a conceptual problem resides in the values 
of the nematic order parameter itself. The thermal ex- 
pansion data sets for all our materials, as well as the 
analogous studies by other groups, e.g. [[l3[ p2| , were all 
found to follow A = 1 + a Q very faithfully, but with Q 
varying in the manner of a virtual second-order phase 
transition, Q oc \T — T ni |^ |2(| . In our present materials 
the critical exponent of the apparent critical behaviour, 
£, takes values between 0.14 and 0.31 (Table [□]). There 
are several possible reasons for this observation of an ap- 
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FIG. 3: Uniaxial thermal expansion of nematic elastomers. 
Plot (a) shows the raw data for the extension of the freely 
suspended elastomers. With plot (b), the temperature for 
each sample is scaled with the extrapolated nematic-isotropic 
transition temperature T n i. Below this point, each data set 
is fitted with a critical power law of the form (A — 1) oc |T — 
T n i|^. The main-chain containing sample Azol8MC shows a 
pronounced pre-transitional expansion (see text). 



parent critical behaviour in a system that, by symmetry 
arguments, should undergo a discontinuous first order 
transition. One likely explanation is that the kinetic re- 
sponse of the rubber-elastic network is such that it is 
unable to perform the discontinuous jump required by 
the changing underlying nematic order parameter. Such 
an explanation, although appealingly simple, has sev- 
eral weaknesses, the main one being that practically no 
dependence on cooling/heating rate is observed, which 
should be expected if the kinetic retardation is at work. 

Another possibility is that in forming the monodomain 
nematic materials, an effective aligning molecular field 
is frozen in by crosslinking. This makes the first-order 
nematic-isotropic transition supercritical and, essentially, 
eliminates the true isotropic phase altogether: even at 
high temperature there would still be a small residual 
order remaining. This could still be the case, even though 
in our materials we could not experimentally detect any 
residual order above T n i and observe a clear and well- 
defined heat signature of the transition on DSC. 

A perhaps more intriguing possibility is the presence 
of quenched sources of random disorder in the nematic 
field (the network crosslinking points), which genuinely 
transforms the weak first order transition into a contin- 
uous one, an effect already discussed in the context of 
spin glasses |27j. The present paper does not attempt 
to address this issue. We focus instead on a different 
physical property of a photo-sensitive nematic elastomer 
and simply use this experimentally observed, continuous 
"critical" variation of Q(T) as the foundation for describ- 
ing the UV-induced response of such an elastomer. 

It seems that a consistent value of the "critical expo- 
nent" is obtained, £ ~ 0.2. The anomalously small value 
of £ for Azol8MC can be attributed to the significant dif- 
ference in this sample's response time from all the others: 
kinetic effects are markedly relevant and the long DiMC 
chains have been shown to slow down, and even com- 
pletely freeze their dynamics, when in the hairpin regime 
of their nematic conformation p8|. This would result 



TABLE II: Experimentally determined transition points, T n i 
(used as a parameter in fitting by the power- law), and the 
apparent critical exponents, £, for all materials studied. The 
samples are ranked by their overall azo-content. 



Sample 


T (0) ( o C) 




DiAzolOO 


120.3 


0.31 


Azo30 


76.8 


0.21 


Azol8 


67 


0.19 


Azol8MC 


75.2 


0.14 


DiAzo9 


66.5 


0.20 



in the non-equilibrated thermal expansion of the sample 
lagging the actual, equilibrated order parameter. This 
sample also shows a very significant pre-transitional ef- 
fect, with the residual nematic order and apparent sam- 
ple extension persisting far above the critical tempera- 
ture. This phenomenon has been examined in a separate 
study [|l0| and is due to a composite nature of this ma- 
terial in which the side-chain and main-chain polymer 
chains are, crudely, in equal proportion. Another point 
of interest is the presence of highly ordered phases in 
DiAzolOO and Azo30 at lower temperatures, where the 
dependence of Q, and hence A, on temperature deviates 
from the attributed nematic behaviour, and becomes es- 
sentially constant. However, this will not be dwelt upon, 
as this work focuses exclusively on the nematic phase and 
the nematic-isotropic phase transition. 

The observed dependence of Q and A on temperature is 
very important in analysing and modelling the response 
of the liquid crystal elastomers to UV-irradiation. In 
doing so, one sample (Azol8), which proved to have the 
most pronounced response to UV-irradiation, is studied 
in greater detail in Section |y[ 



B. Response to UV-irradiation 

In order to compare different compositions and 
crosslinking topologies of the nematic elastomers in their 
response to UV-irradiation, all samples were equilibrated 
before irradiation at approximately 10°C below their crit- 
ical temp eratur e T ni , determined by the thermal analysis 
in section III A. Such proximity to T n i also aims to max- 



imise the variation of A, whilst avoiding the singularity 
at the phase transition. 

During irradiation, the sample temperature was, with- 
out exception, found to increase in an approximately ex- 
ponential manner, mirroring the behaviour of the sponta- 
neous uniaxial extension, A. Similarly, as the elastomers 
relaxed after the radiation was switched off, the mea- 
sured temperature again mirrored A, decreasing back to 
its pre-irradiated value. In both cases, this exponential 
variation in temperature was found to be modulated by a 
small-amplitude oscillatory behaviour. Possible explana- 
tions for this thermal phenomenon involve a trivial heat- 
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FIG. 4: The response of Azol8 to UV irradiation. Both the 
relative length of the sample (uniaxial strain A, measured 
with respect to the isotropic state: filled symbols) and the 
measured temperature (open circles) are shown. The ±0.1° C 
errors in temperature have been omitted for clarity. Solid 
lines are fits of both the on and off branches by a simple ex- 
ponential, which at this point should be regarded as only a 
guide to the eye. 



ing effect from the UV lamp interfering with the active 
temperature control of the oven. However, such an ad- 
ditional heating power should not result in the overall 
rise of saturated temperature value, as seen in Fig. ||: 
the oven would be expected to reduce its own heating 
and bring the temperature back to its set value. Also, it 
would be an inconceivable coincidence if the characteris- 
tic rate of the oven heating system are exactly the same 
as that of all our different materials. 

By closely examining the Fig. || it becomes clear that 
we observe the effect of spontaneous sample heating, oc- 
curring due to the exothermic nature of isomerisation of 
azobenzene groups. In a following section we shall de- 
scribe the population dynamics between the trans and 
cis N=N isomers. The relative proportion of these states 
is controlled by the rates of forward and reverse reac- 
tions, which reach the balance in the saturated state (af- 
ter a long time of irradiation). If one assumes that the 
heat release on the return to the trans ground state is 
proportional to the rate of change of cis concentration, 
h C i S , then the observation in Fig. [| becomes sensible. Ini- 
tially, all azobenzene groups are in their trans state and 
the sample is thermally equilibrated with the oven (at 
59.5°C in this case). On irradiation, the increasing num- 
ber of cis isomers is created, and also accompanied by 
an ever increasing relaxation flux back to the trans state. 
This heats the sample, at the rate exactly matching the 
rate of isomerisation itself. Finally, a saturation state is 
reached, when the UV-stimulated trans^cis flux is equal 
to the relaxation flux: at this stage the heat release in 
the sample is constant and maximal. On switching the 
UV-light off the exothermic cis^ trans relaxation contin- 
ues to heat the sample, but with a decreasing power as 
the overall concentration n C i S decreases. 

All the samples responded to UV-irradiation in a man- 
ner similar to Azol8, illustrated in Fig. [|, differing only in 



TABLE III: Sample contractions upon UV-irradiation, rela- 
tive to the expected thermal contraction due to the accom- 
panied temperature increase above the starting point (last 
column). The increase of this ration above unity is a genuine 
measure of the effect of UV isomerisation. 



Sample 


AA/AA(T) 


^exp ( C) 


DiAzolOO 


1.21 


110.5 


Azo30 


1.49 


66.9 


Azol8 


2.82 


57 


Azol8MC 


1.03 


65.0 


DiAzo9 


1.63 


56.8 



the magnitude and speed of their response. The resultant 
contraction upon irradiation for each sample, relative to 
the expected contraction due to the aforementioned in- 
crease in temperature, are summarised in Table III. As 
can be seen, the proportion of UV-isomerisable azoben- 
zene rods in the sample is by no means a direct indication 
of its response to UV-irradiation, interpreted as the max- 
imal sample extension achieved on irradiation. 

Azol8 gave the largest response to radiation, but 
the very similar Azo30, with its higher proportion of 
azo-rods, gave a significantly smaller response. This 
must be due to the position of the azo-group within 
the corresponding mesogenic moiety - with Azo30, it 
is only weakly coupled to the polymer backbone by a 
long spacer, in contrast to Azol8, where a strong par- 
allel alignment between the rod and the backbone is in- 
duced. This could result in the kinked cis-form of the 
moiety having less of a destabilising effect on the back- 
bone anisotropy, and so the elastomer response to UV- 
irradiation is curtailed. 

DiAzo9 is very similar to the material studied by 
Finkelmann et al. [15 . Despite containing less than 
a third of the proportion of azobenzene rods found in 
Azo30, its response to UV-irradiation is comparable in 
magnitude, which should be attributed to the crucial po- 
sition of the photo-isomerisable group, in the centre of 
the predominant, short crosslinker. 

Azol8MC might have been expected to have a very 
large response compared to the other materials due to 
its large backbone anisotropy and large coupling con- 
stant to the underlying order parameter, illustrated by 
its large thermal expansion amplitude. However, this 
was definitely not the case for UV irradiation. One may 
speculate the cause may be not only that the actual pro- 
portion of azobenzene rods was quite low for this sample 
(much of the volume is occupied by the main-chain poly- 
mer crosslinker), resulting in a smaller absolute change 
in the order parameter, but more importantly, that large 
DiMC chains can essentially freeze the elastomer dynam- 
ics in response to changes in Q [ p8| , quenching its mani- 
festation through A. 

All the mesogenic groups and crosslinkers in DiAzolOO 
were photo-sensitive, and so one might have expected 
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an exceptional mechanical response to UV-ir radiation. 
However, the critical temperature, T n j, of this sample 
was much higher than that of the other materials, and 
to maintain consistency between experiments on differ- 
ent samples this required an initial equilibrium temper- 
ature of approximately 110°C. At such a temperature 
the thermally-activated relaxation taking the azobenzene 
groups from their metastable kinked cis-shape back to 
the rod-like trans-shape would be much accelerated. As 
a result, one is perhaps unable to isomerise a signifi- 
cant proportion of azobenzene rods at any given moment 
of time. A lower equilibrium population of cis-isomers 
would then result in a smaller overall change in Q, and 
hence in A. 

The solid curves in Fig. [I] are simple exponential fits 
to the data, for both the UV-on and the UV-off states. 
Apparently, it is a reasonably good fit. However, the 
parameters are rather arbitrary and appear to be quite 
unrelated to the corresponding spontaneous thermal ex- 
pansion curves. It will shortly be seen that a different, 
universal theoretical dependence on temperature and ir- 
radiation time describes the data with a much greater 
degree of fidelity. Accordingly, the solid lines in Fig. || 
serve only as a rough guide to the eye. 



IV. MANIPULATION OF NEMATIC ORDER 

A. Population dynamics 

The internal energy of rod-like trans, and the kinked 
cis mesogenic moieties is adequately described as a two- 
level potential well, with the former configuration ener- 
getically more stable than the latter (Fig. ||). In the 
experiments conducted, there are three possible ways of 
moving between states, each characterised by a separate 
time constant: 

• Thermal excitations up, from trans to cis, with the 
characteristic time given by r tc = roe u ^ kT , with To the 
bare attempt rate. 

• Thermally driven relaxation down, from cis to trans, 
with the time constant r ct = roe A / fcT . 

• UV-isomerisation, occurring with the rate l/r uv = V oc 
radiation intensity. 

One expects that the barrier height U is much greater 
than the metastable state depth A, and hence one can 
neglect the rate of spontaneous, thermally driven trans 
to cis transitions, although certainly not its reverse. 

Let us assign the total (fixed) number density of photo- 
isomerisable rod-like groups in a given elastomer as no = 
ritrans + n C i S , with the only independent variable defined 



as n 



n tr 



so that n r 



no — n. Thus, the rate of 



change in the population of trans isomers can simply be 
written as; 

d 11 

—n(t,T) = -r\n n H (n - n) (1) 

ot T tc r ct 

In the above approximation, the relaxation time r tc is as- 




FIG. 5: Schematic depiction of the potential energy profile 
of a photo-sensitive moiety verses configuration, showing the 
relative stabilities of the cis and trans isomers. 



sumed to be very large, and can therefore be neglected. 
The appropriate initial condition is then n t =o = (that 
is, all azobenzene groups are in their rod-like trans con- 
formation) and the dynamic solution to Eq. |l| is simply; 

n{t, T) = n ° [1 + T ctV exp (-f/reff)] (2) 

J- + T c tT] 



with T eff = (77 + l/r ct ) 



Hence, one can see that the increasing concentration of 
kinked cis-isomers, which act as impurities for the ne- 
matic ordering, will eventually reach a plateau due to 
the opposing thermal back-reaction, itself characterised 
by the rate l/r ct . Both the rate and the final saturation 
value of the isomer density are crucially dependent on 
the irradiation rate 77; this could be different for different 
materials, sample shapes and experimental conditions as 
it is affected, among other factors, by the radiation ab- 
sorption in the medium. 

The relaxation process, after switching the radiation 
off (that is, setting 77 = 0), is the thermally driven flux of 
ds-isomers into the trans ground state over the barrier A. 
Assuming a fully saturated irradiated state, with n C i S = 
no — — r C f 77/(1 + r ct r]) at t = (and still neglecting 
the spontaneous thermal trans^cis process), results in 
the relaxation law 



n(t,T) = n 



r c tV 

1 + T ct r] 



exp (-t/r ct ) 



(3) 



Naturally, the relaxation time scale is purely r c t, a prop- 
erty of an individual azobenzene group, which should re- 
main at least approximately universal between all differ- 
ent samples under study. 



B. Macroscopic shape and nematic order 

The influence of impurities on the nematic order is well 
studied [^9|. One of the simplest ways of representing 
their effect is via the Landau-de Gennes mean-field for- 
mulation, which describes the local free energy density 
as a power series in the magnitude of the scalar order 
parameter, Q. Such an approach is commonly used in 
the thermodynamics of second order phase transitions, 
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0.2 0.4 0.6 0.8 

Reduced temperature T/T ni 

FIG. 6: The theoretically calculated effect of increasing the 
intensity of UV-irradiation (measured by the effective non- 
dimensional rate rof]) on the uniaxial thermal extension of a 
photo-sensitive liquid crystal elastomer, according to Eq. 0. 



as well as in application to liquid crystals. Although it is 
much less suitable for the nematic-isotropic transition in 
liquid crystals, which should be first-order by symmetry, 
many useful qualitative results are obtained within the 
Landau-de Gennes formalism. 

However, in view of many limitations of such a phe- 
nomenological mean-fi eld t heory, we adopt a different 
approach. In section [II A, the actual form of the de- 



pendence of the order parameter on temperature was es- 
tablished experimentally as being a virtual, continuous 
transition with the apparent critical exponent £ quite 
far from the mean-field value predicted by the Landau- 
de Gennes approach (and certainly appears a contin- 
uous transition). Therefore, we simply take the em- 
pirical dependence, Q ex \T — T n i|^, and modify it by 
substituting in the impurity-shifted critical temperature, 
T ni = T ni [n(t,T)]. 

The effect of impurities on the critical temperature is, 
in the first approximation, a linear function of the pop- 
ulation of ds-isomers. The size of this population, n C i Sl 
defined above as no — n(i), is itself a function of time; 



T ni (n) 



where T c = is the transition temperature in the 

"pure" material containing only the trans rod-like units 
(given in the Table || for our materials). Combining this 
with the known dependence of the spontaneous uniaxial 
sample extension, A, on Q, one obtains the following re- 
lation to describe the observed mechanical response to 
UV-irradiation 



(3[n - n(t)], 



(4) 



A U v(t,T) 



= 1 + a 



(5) 



Pn 1 



TctV 

1 + r ct rj 



(1 



-*/ 



T e ff \ 



T 



Plotting A against temperature for different values of the 
intensity parameter 77, Fig. ^ one sees that the expected 
consequence of UV-irradiation is simply to lower the criti- 
cal temperature of the nematic-isotropic phase transition. 
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Time (s) 

FIG. 7: The response of Azol8 to UV-irradiation, and the 
corresponding theoretical fit by Eq. ([]), solid line. Note that 
even the minor intrinsic temperature oscillations are, in fact, 
directly manifested in A(£,T), and are almost perfectly de- 
scribed by the theoretical model (see text). 



The relaxation of sample shape after switching the ra- 
diation off is determined by the rate of recovery of the 
initial population of trans isomers, according to Eq. (||); 



A off (t,T) = l 



f3n 



r c tV 

( 

1 + r ct rj 



-t/r c 



-it 



, (6) 



where, as in the Eq. (|), we have assumed that a fully 
saturated state was achieved on preceding irradiation. 

One important parameter has to be used in the further 
analysis: the relaxation time scale of the cis^ trans pro- 
cess. From the data in Fig. || for the Azol8 at ~ 59°C it 
is estimated as r ct « 1479 s, assuming that the dynamics 
is dominated by the rate of isomerisation (and not, for in- 
stance, the viscoelastic relaxation of a rubbery network). 
In the following analysis of experimental data, the effect 
of (Arrhenius) temperature variation of r ct within a given 
data set are, in the first approximation, neglected. This 
is justified by the relatively narrow range of absolute tem- 
peratures over which a given data set varies. However, 
across data sets, r ct is shown to vary significantly with 
temperature. 



C. Nematic destabilisation by UV radiation 

Applying the Eq. (||) to the data in Fig. |], one obtains 
a remarkably improved fit (Fig. Q). In analysing the data 
we take the fixed values for the coefficient a = 0.22 and 
the exponent £ = 0.19, as deter mined by the thermal ex- 
pansion measurements of section III A and the relaxation 
time T ct = 1479 s, as estimated from a crude exponential 
fit of Fig. ||. The values of current temperature T are 
taken from the corresponding experimental data points. 
The free parameters are, thus, only the radiation rate 



r] => 2.2 • 10 s and the product (3 no 11.8, which is 
a measure of the overall concentration, and the efficiency, 
of azobenzene groups. 
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The result of such a fit for A(t, T) is very satisfying. 
The initial quasi-linear contraction and its subsequent 
rapid flattening and saturation are fully described, unlike 
the rough exponential fit in Fig. [|. More importantly, we 
now discover that what looked like noise in both the tem- 
perature and the strain \(t) in the initial graph, Fig. U, is 
in fact a consistent small-amplitude temperature oscilla- 
tion (probably due to the active oven control operation), 
which is also fully reproduced in the mechanical shape of 
the sample. This is a rather stringent test of the theoreti- 
cal model and the parameters used to fit the experimental 
data. 

All the other data sets for Azol8 at different temper- 
atures are fitted equally well, as shown in Fig. ||. Once 
again, we fix the value of strain-order coupling coeffi- 
cient a = 0.22 and the exponent £ = 0.19, as deter- 
mined by the thermal expansion measurement for this 
sample, and the relaxation time r ct = 1479 s. On phys- 
ical grounds, we also insist that the coefficient /3no re- 
mains constant (~ 12) between these data sets, since it 
is determined by the intrinsic properties of the nematic 
order and the given concentration of azobenzene groups. 
With these stringent constraints, one obtains the values 
of the only remaining free parameter, the effective rate of 
induced isomerisation. The inverse of 77, the characteris- 
tic time for the photo-induced transacts isomerisation, 
r uv = I/77, is plotted in Fig. |)](a) and has an average 
value r uv ~ 3 • 10 3 s. Although there is a certain spread in 
these values, we find the outcome of this analysis quite 
satisfactory given the number of constraints imposed on 
the fitting, certain error in the data and, in spite of that, 
a clear evidence that not only the overall trend but also 
the small oscillations of temperature are well reproduced 
bytheEq. (|). 

The values of the characteristic time for the reverse 
thermal cis^ trans relaxation, r c t, are consistent with 
those reported in the literature for azobenzene photo- 
chemistry |p5| , as well as within themselves, Fig. ||(b). 
The crude estimate of the cis^ trans energy barrier from 
the Arrhenius fit of the relaxation time returns a reason- 
able value A - 4 • 10" 20 J - 10/cT. This is in contrast 
with the isomerisation time r uv , which was found to be 
a constant, or very slowly decreasing function of temper- 
ature. However, the errors in determining n were large, 
and there is much room for further studying its depen- 
dence on temperature. 



D. The vicinity of critical point 



transition line into the isotropic phase. Continuing the 
study of Azol8, we now subject it to UV-radiation at a 
temperature T exp ~ 63° C, only 4° below T c . 

Fig. [l(] indicates that the mechanical response of the 
sample on irradiation, A(t), starts exactly in the same 
way as at all other temperatures. However, the initial 
steep drop in A (and in the underlying nematic order 
Q) rapidly reaches the point of A = 1 (or Q = 0); the 
material becomes fully isotropic and no longer has any 
macroscopic response to the further increase of the popu- 
lation of cis-isomers. The irradiation continues for some 
time, as illustrated by the continuing rise in tempera- 
ture, which nevertheless has never reached the "bare" 
T c . Then the UV-light is switched off and the tempera- 
ture begins to fall, in the same way as all previous mea- 
surements have shown. However, nematic order is not 
established for some time: clearly the concentration of 
czs-impurities has to decrease sufficiently to allow the ne- 
matic phase to appear at a given temperature. 

Accordingly, there are two simultaneous thermal ef- 
fects in the material. One is the rise (in the UV-on state), 
or fall (in the off state), of the actual temperature of the 
sample. The other is the underlying shift of the critical 
temperature of the nematic phase transition, described 
by the Eq. (||) with the time-dependent concentration 
n(t) given by Eqs. (|^) and (|^), respectively (in the latter 
case of thermal czs— » trans relaxation, Eqs. (||) and @ 
assume that the state of full saturation has been reached 
on irradiation, otherwise a different initial condition has 
to be used, altering the prefactor to the exponential). 
It is straightforward to apply the corresponding analysis 
and determine the time when, in the UV-on state, the 
decreasing critical temperature T n i = T c — /3[no — n(t)] 
becomes equal to the increasing sample temperature. 
Equally, in the W-off state, we find the time needed 
for the increasing T n i to intersect the decreasing actual 
temperature T. 

Applying this analysis to Azol8 in Fig. |l0|, one ob- 
tains that the nematic order disappears at t « 1100 s and 
subsequently re-enters the nematic phase, is t « 4850 s. 
Both times are marked with arrows in Fig. [l0[ clearly 
the agreement with the observed change in the sample 
length, A = L/Lq, is excellent. This is another confir- 
mation of the theoretical model, which assumes simple 
thermal activation laws for characteristic times, the lin- 
ear dependence of critical temperature on the concentra- 
tion of ds-isomers and, much less trivially, the empiri- 
cal critical dependence of the nematic order parameter 
Qoc \T-T*\Z. 



Since the effect of UV radiation on a nematic rub- 
ber containing isomerisable azobenzene groups is clearly 
identified as an effective shift of critical temperature of 
the underlying nematic-isotropic phase transition, the 
question arises - what would be the response when the 
material is irradiated at a temperature close to the "bare" 
critical point T c ? In this case, the increasing concentra- 
tion of ds-impurities should take the system past the 



V. CONCLUSIONS 



The derivation of Eq. (ph in section EVB, cumulat- 



ing in the analysis of relaxation data sets in Figs. |7|-|g, 
makes it clear that, once the number of kinked moieties 
reaches a steady-state population, the resultant effect of 
UV-irradiation is equivalent to simply reducing the criti- 
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FIG. 8: The response of Azol8 to UV-irradiation at different starting temperatures, ~ 57, 50 and 25°C for plots (a), (b) and 
(c), respectively. The corresponding description of A(t,T) with the theoretical model, Eq. (Q), is shown by the solid lines. The 
relative error increases for lower temperatures due to the smaller overall magnitude of the UV-response. 
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FIG. 9: (a) The irradiation time constants, I/77, at different 
temperatures, are all of the same order of magnitude and show 
a constant or slowly decreasing dependence on temperature. 
The circled point, near T n i, h as an ambiguously high error 
associated with it (see section IV D). (b) Relaxation times 
r c t gleaned from a simple exponential fit of the relaxation 
data. The activation temperature dependence allows a crude 
estimate of the isomerisation energy barrier A~4-10 _20 J~ 
lOfcT. 




1.0 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
20 30 40 50 60 70 80 90 

Temperature (°C) 

FIG. 11: The asymptotic values of A(t,T), for Azol8 irradi- 
ated for a long time at several fixed temperature, after cor- 
recting for the temperature increase associated with irradia- 
tion, f3[no — n], which correspond well to the idea of a shifted 
critical temperature T*, as illustrated theoretically in Fig. n 



1.25-^UV On- 



rUV Off 




4000 6000 
Time (s) 



10000 



FIG. 10: The response of Azol8 to UV-irradiation in the 
vicinity of critical point. Arrows indicate the moments of 
time when the material first becomes fully isotropic on irra- 
diation, and then re-enters the nematic phase on relaxation 
(the associated measured temperature of ~ 64° C is clearly 
insufficient to account for this transition). 



cal temperature of the nematic-isotropic phase transition 
by an amount 



AT = P[n -n oo (T)}=Pn 



TctV 
1 +T ct T) 



which is dependent on the intensity of the ultraviolet 
light used (through 77). Hence, a slightly modified form 
of A = 1 + aQ, making only the replacement T n [ — > 
T*, should fully describe the final, steady-state length of 
the sample, X uv (t — > 00), as a function of temperature. 
All the other parameters in Eq. (||), in particular the 
exponent £, should remain unchanged. Indeed, one finds 
that this idea does describe X uv (t — > 00) very well, as 
shown in Fig. |ll|. The corresponding value of shift AT, 
which best fits the data in Fig. [□], is found to be (3.3 ± 
0.7)°C which agrees well with the average asymptotic 
value of /3[no — n^] gleaned from extrapolating the data 
sets in Figs. ^1 of (3.5±1.9)°C. 

This work, has studied a wide range of photo- 
sensitive materials, based on azobenzene molecular moi- 
eties, whose differing compositions and topologies are 
reflected in their varied mechanical response to irradi- 
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ation with ultraviolet light. The fundamental details 
of this response were investigated, in particular its ki- 
netics and the variation of the critical temperature near 
the nematic-isotropic phase transition. The obtained ex- 
perimental results are very well described by a theoret- 
ical model assuming a simple linear shift in the critical 
temperature with the concentration of kinked cis groups, 
which results in a rather non-trivial model for the depen- 
dence of the underlying nematic order parameter Q and 
macroscopic uniaxial mechanical extension A on irradi- 
ation time and temperature. Basic material parameters 
and constants were also determined, which will allow sub- 
sequent investigations to further probe the details of this 
fascinating phenomenon. With such a novel, non-invasive 



method of control over mechanical actuation, there is a 
great potential for application, but still very much is left 
unknown. 
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